Size Dependence and Convergence of the Retrieval Parameters of Metamaterials 
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We study the dependence of the retrieval parameters, such as the electric permittivity, e, the 
magnetic permeability, fi, and the index of refraction, n, on the size of the unit cell of a periodic 
metamaterial. The convergence of the retrieved parameters on the number of the unit cells is also 
examined. We have concentrated our studies on the so-called fishnet structure, which is the most 
promising design to obtain negative n at optical wavelengths. We find that as the size of the unit 
cell decreases, the magnitude of the retrieved effective parameters increases. The convergence of 
the effective parameters of the fishnet as the number of the unit cells increases is demonstrated but 
found to be slower than for regular split ring resonators and wires structures. This is due to a much 
stronger coupling between the different unit cells in the fishnet structure. 
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INTRODUCTION 



The recent development of metamaterials with nega- 
tive refractive index has confirmed that structures can 
be fabricated that can be interpreted as having both a 
negative effective permittivity, e, and a negative effec- 
tive permeability, fi, simultaneously. Since the original 
microwave experiments for the demonstration of nega- 
tive index behavior in the split ring resonators (SRRs) 
and wires structures, new designs have been introduced, 
such as the fishnet, that have pushed the existence of the 
negative refractive index at optical wavelengths P, [3, [1] ■ 
Most of the experiments [3, 0, 0, S| with the fishnet 
structure measure the transmission, T, and the reflection, 
R, and use the retrieval procedure 0, [3, U, 13, l3| to 
numerically obtain the effective parameters, e, and n. 
The only direct experimental way that was able to mea- 
sure both the phase and the group velocity, was done by 
pulse measurements [6] . In all the other experiments only 
T and/or R are measured and then the retrieval proce- 
dure is used to obtain the effective parameters. There 
is an unresolved issue in the retrieval procedure, what is 
the size of the unit cell that needs to be used. Most of 
the retrieval procedures use only one unit cell and made 
the assumption that if more than one unit cell will be 
used, the results will be the same as for the single unit 
cell. This is an untested assumption and Koschny et al. 
[11] have systematically studied this assumption for the 
SRRs and wires. It was found that this assumption is 
correct if the wavelength is 30 times larger than the size 
of the unit cell. No systematic study have been done on 
the dependence of the retrieval parameters on the size 
of the unit cell of the fishnet structure. In this work, 
we tackle this problem and also study the convergence of 
the effective parameters from the retrieval procedure as 
the number of the unit cells increases. We also, exam- 



ine how the figure of merit (FOM), which is defined as 
— Re(n)/Im(n), changes with the size of the unit cell and 
the number of the unit cells. 



NUMERICAL CALCULATIONS 
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FIG. f : Schematic of a fishnet structure (a); a single unit cell 
with geometric parameters marked on it (b). 

Our numerical simulations were done with CST 
Microwave Studio (Computer Simulation Technology 
GmbH, Darmstadt, Germany), which uses a finite- 
integration technique, and Comsol Multiphysics, which 
uses a frequency domain finite element method. The 
schematics of the unit cell of fishnet structure are shown 
in Fig. Notice that he propagation direction is 

perpendicular to the plane of the fishnet. In Fig. [5] we 
present the results for the transmission, T, reflection, R, 
and absorption, A — 1 — T — i?, for a given size of the 
unit cell. 

In Fig. [31 we present the results of the real part of the 
retrieved effective parameters, n, fj, and e, as a function 
of the frequency. Notice that as the size of the unit cell is 
getting larger the strength of the resonance of the mag- 
netic permeability fi is getting weaker. The same is true 
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FIG. 2: Simulated transmission (red solid), reflection (blue 
dashed) and loss (green dotted), which is proportional to ab- 
sorption, of a fishnet structure. The geometric parameters 
are = ay = 15 mm, = 4 mm, Wx = 6 mm, Wy = A mm, 
s = 0.5 mm and t = 0.05 mm and the dielectric constant of 
the spacer is e = 5. 



for the magnitude of the index of refraction, n, which is 
also negative. For the dimensions we have chosen, the 
resonance frequency is 9.9 GHz, which is equivalent to 
a wavelength of A = 30.3 mm. So the ratio of A/a^ ~ 
30, 15, 7.6 for the three sizes {a^ =1 mm, 2 mm and 4 
mm) of the unit cell examined in Fig. [3] The real part 
of e is always negative, and shows an anti-resonance [loj 
corresponding to the magnetic resonance at 9.9 GHz. 

Another important issue [l^ that the community of 
metamaterials has not completely understood is why the 
width of the negative n is wider than the width of the 
negative ^, as can be seen in Fig. [3]Ja) and[3fb). In Fig. 
m we plot together the fx and n (real parts) around the 
magnetic resonance frequency of 9.9 GHz for the = 
1mm case. Notice that /x < for frequencies higher than 
9.9 GHz. However n is also negative below 9.9 GHz where 
fi > 0. So we have two regions that give negative n. The 
one above 9.9 GHz (drawn as blue in Fig. |4]) which has 
both e and /x negative and therefore n is also negative. 
The other region, below 9.9 GHz (shown as yellow in 
Fig. m has e < but ^ > and still give n < 0. The 
reason is that one can obtain Re(n) negative provided 
that Im(e)Im(z) > Re(e)Re(2;), where z = \/u/e. For 
the details of this explanation see Ref. [il,|iB,[li. The 
normalized loss, (1 — T — i?)/(l — R), which is the ratio 
between the absolute loss, 1 — T—R, and the total energy 
entering the fishnet structure, 1 — i?, is also shown in 
Fig. m The normalized loss is a better measure of the 
absorption in a medium, since it removes the effect of 
the reflection due to the impendence mismatch at the 



interface between the medium and the vacuum space, 
which results in lower absolute loss. If we compare the 
frequencies with same Re(n) value in the regions below 
9.9 GHz (yellow) and above 9.9 GHz (blue), one can see 
that the normalized loss in the region below 9.9 GHz 
(yellow) is much higher than the region above 9.9 GHz 
(blue). The transmission, T, as shown in Fig. [21 in region 
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FIG. 3: Retrieved real part of refractive index, n, (a), perme- 
ability, fj,, (b), and permittivity, e, (c), from simulated data 
shown in Fig. [J] using unit cell size in propagation direction 
az = 1 mm (red solid), =2 mm (blue dashed) and = 4 
mm (green dotted). 



above 9.9 GHz (blue) is much higher than in the region 
below 9.9 GHz (yellow). Therefore, the frequency region 
above 9.9 GHz (blue) is the frequencies range which is 
suitable for the applications of the negative refractive 
index medium. 

In Fig. [5l we present the results of the figure of merit 
(FOM=— Re(n)/Im(n)) as the size of the unit cell of the 
fishnet structure increases. Notice that as the size of 
the unit cell is getting smaller, the FOM increases and 
reaches a value of 6.5, when the size of the unit cell along 
the propagation direction is equal to 1 mm. Notice also 
that the maximum of the figure of the merit is higher than 
the resonance frequency, which is 9.9 GHz. The reason 
is that the Im(n) takes its larger value at the resonance 
frequency and decreases as one goes away from the res- 
onance frequency. This is the reason that the FOM=6.5 
for az — 1 mm has its maximum value at 9.95 GHz. 

In Fig. [6l we present the results of the retrieved n, 
/i and e (real parts only), as a function of frequency, 
for THz frequencies. As we have discussed above, the 
fishnet structure is the only structure that has been 
proven to give negative n at THz and optical frequen- 
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FIG. 4: Retrieved real part of effective refractive index, Re(n), 
(red solid), permeability, Re(/i), (blue dashed), and normal- 
ized loss, {1 — T — R)/{1 — R), (green dotted) (normalized 
loss has been magnified by a factor of 5 to improve visibility) . 
The shadow regions, below 9.9 GHz (yellow) and above 9.9 
GHz (blue), show frequency regions with Re(n) < —0.5. 
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5: Figure of merit, FOM 
= 1 mm (red solid) 
— 4 mm (green dotted). 



— Re(n)/Im(n), for unit cell 
— 2 mm (blue dashed) and 
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FIG. 6: Retrieved real part of effective refractive index, n, (a), 
permeability, /i, (b) and permittivity, e, (c) from simulation 
data of a fishnet structure with parameters = 200 nm (red 
solid), az = 300 nm (blue dashed), and = 400 nm (green 
dotted). The other geometry parameters are = ay = 600 
nm, Wx ~ 400 nm, Wy = 200 nm, s = 100 nm and t = 60 nm, 
and the dielectric constant of the spacer is e = 1.9. 



cies 0, 0, [3 ■ To obtain negative /i and negative n at 
THz frequencies, we have to scale down the dimensions 
given in Fig. [51 The size of the unit cell along the propa- 
gation direction — 200 nm, 300 nm and 400 nm. The 
resonance frequency is approximately 235 THz, which is 
equivalent to a wavelength A — 1.28 /im. So the ratio of 
A/flz =6.4, 4.27 and 3.2 for the three sizes Qz =200 nm, 
300 nm, and 400 nm of the unit cell examined in Fig. 
[6l Notice in Fig. [H^a) that as the size of the unit cell 
increases the retrieved n goes from negative to positive 
values. The reason for the positive value of n is the fol- 
lowing: As the size of the unit cell along the propagation 
direction increases, the effective plasma frequency of the 
fishnet structure, as shown in Fig. El^c), decreases and 
becomes smaller than the magnetic resonance frequency 
of fi, as shown in Fig. El^b). Therefore, we don't have an 
overlap of the negative region of e and fi. 

In Fig. [71 we present the results of the transmission 
of the one unit cell, four unit cells, eight unit cells and 
ten unit cells (along propagation direction) for a fishnet 



structure that has a magnetic resonance at 15.65 GHz. 
Notice that as the number of the unit cells increases, the 
transmission results converge to a common value. This 
is seen more clearly in Fig [8l where we present the re- 
trieved results for the effective index of refraction. One 
can clearly see that the results for eight and ten unit cells 
are exactly the same. This convergent length dependence 
is a required property to attribute effective medium be- 
havior to a metamaterial ll| . Notice that the converged 
results for the effective parameters are different from the 
effective parameters results of the one unit cell. As one 
can see from Fig. \8la) , the retrieved results for the effec- 
tive parameter n for more than one unit cells are com- 
pletely different for the results of the one unit cell. This 
is true, especially when the Re(n) is relative large. At 
the high frequency edges of the negative n region, there 
is no much difference between the results of the one unit 
cell with the results of more than one unit cells. In Fig. 
[Sl[b), one notices that the FOM does not change dra- 
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FIG. 7: Simulated transmission spectra of one layer (red 
solid), four layers (blue dashed), eight layers (green dotted) 
and ten layers (black dash-dotted) of the fishnet structure. 
The geometric parameters are Ux = ay = 10 mm, = 4 mm, 
Wx = 4: mm, Wy — 3 mm, s — 1 mm and t = 0.05 mm, and 
the dielectric constant of the spacer is e = 4. 



matically as one uses more unit cells to do the retrieval 
procedure. As one uses more unit cells, the FOM has 
more oscillations as a function of frequency. 
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FIG. 8: Retrieved real part of effective refractive index, n, and 
figure of merit of the simulated structure described in Fig. [7] 



Fig. [9] shows the real parts of effective refractive in- 
dex, Re(n), as a function of frequency for one layer and 
two layers fishnet structure with different unit cell size: 
az ~ 1 mm, = 2 mm and Cz = 4 mm. The reso- 
nance frequencies, i.e. frequency with maximum |Re(n)|, 
are 9.926 GHz (red dashed), 9.896 GHz (blue dashed) 
and 9.894 GHz (green dashed) for one layer; 9.973 GHz 
(red sohd), 9.927 GHz (blue solid) and 9.900 GHz (green 
solid) for two layers. The difference in value of resonance 
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FIG. 9: Retrieved real part of refractive index, n, from simu- 
lated data using unit cell size in propagation direction = 1 
mm (red), = 2 mm (blue) and = 4 mm (green). Both 
one layer (dashed) and two layers (solid) results are shown. 



frequencies between one layer and two layer simulations 
are 77 MHz (red, with = 1 mm), 31 MHz (blue, with 
az — 2 mm) and 6 MHz (green, with = 4 mm). It 
is shown that the resonance frequency, /m, of two lay- 
ers results decreases as increases, while fm of the one 
layer only shifts slightly, when the size of the unit cell 
increases. The frequency shift of two layers simulation is 
due to the coupling between the two neighboring layers of 
the fishnet. Notice that as the unit cell size increases, 
the coupling between neighboring layers becomes weaker, 
and therefor the effective refractive index, Re(n), of two 
layers approaches the one layer simulation results. 



CONCLUSIONS 

We have systematically studied the size dependence 
of the retrieved parameters, /x, e and n, of the fishnet 
metamaterial structures as a function of the size of the 
unit cell. We found that the retrieval parameters have 
a stronger resonance behavior as the size of the unit 
cell decreases. We have also studied the convergence of 
the retrieved parameters as the number of the unit cells 
increases. We found that the convergence depends on 
the ratio of A/a, where A is the wavelength and a is the 
size of the unit cell. The larger A/a is, the easier the 
convergence. The convergence of the fishnet structure is 
much slower than that observed in SRRs and/or wires 
systems. Finally, we gave an explanation of why the fre- 
quency width where n is negative is wider than the width 
where /i is negative. 
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